Incoherent Neutrinoproduction of Photons and Pions in a Chiral 

Effective Field Theory for Nuclei 

Xilin Zhang* and Brian D. SerolQ 

Department of Physics and Center for Exploration of Energy and Matter 
Indiana University, Bloomington, IN 4^405 
(Dated: March 8, 2013) 

Abstract 

We study the incoherent neutrinoproduction of photons and pions with neutrino energy E v ^ 
0.5 GeV. These processes are relevant to the background analysis in neutrino-oscillation experi- 
ments [for example, MiniBooNE; A. A. Aquilar-Arevalo et al. (MiniBooNE Collaboration), Bhys. 
Rev. Lett. 100, 032301 (2008)]. The calculations are carried out using a Lorentz-covariant effective 
field theory (EFT), which contains nucleons, pions, the Delta (1232) (A), isoscalar scalar (<r) and 
vector (w) fields, and isovector vector (p) fields, and has SU(2) L SU(2) R chiral symmetry realized 
nonlinearly. The contributions of one-body currents are studied in the local Fermi gas approxima- 
tion. The current form factors are generated by meson dominance in the EFT Lagrangian. The 
conservation of the vector current and the partial conservation of the axial current are satisfied 
automatically, which is crucial for photon production. The A dynamics in nuclei, as a key com- 
ponent in the study, is explored. Introduced A-meson couplings explain the A spin-orbit coupling 
in nuclei, and this leads to interesting constraints on the theory. Meanwhile a phenomenological 
approach is applied to parametrize the A width. To benchmark our approximations, we calculate 
the differential cross sections for quasi-elastic scattering and incoherent electroproduction of pions 
without a final state interaction (FSI). The FSI can be ignored for photon production. 
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I. INTRODUCTION 



This paper is a continuing weak offlQ, focusing on neutrinoproduction of photon, and 
pions from nuclei with neutrino energy E v ^ 0.5 GeV. In Refs. 12J, we introduced the 



A resonance as a manifest degrees_o 
quantum hadrodynamics or Qt 
results are discussed in Refs 



freedom to the effective field theory (EFT), known as 



: D jij- 10 1 . (The motivation for this EFT and some calculated 
UH20J.) To calibrate the reaction mechanism on the 
nucleon level, we studied the productions from nucleons The calculations are motivated 
by the fact that the neutrinoproductions of tt° and photons from nuclei (and nucleons) 



are potential backgrounds in neutrino-oscillation experiments (e.g., MiniBooNE [2lH23j). 
Currently, it is still a question whether the neutral current (NC) photon production might 
explain the excess events seen at low reconstructed neutrino energies, which the MicroBooNE 



experiment plans to answer [24( . Moreover, the authors of Refs. [25|-|28j point out the possible 
role of anomalous interaction vertices involving u(p), Z, and the photon in NC photon 
production. So it is necessary to calculate the cross sections for these processes. Here 
by using the QHD EFT, we study incoherent production, in which the nucleus is excited. 
Coherent production with the nucleus being intact is a topic of future work. 1 We will discuss 
the power-counting 2 of the calculations through which we will show that the contributions 
of the anomalous interactions are small in the incoherent NC production of photons (where 
they contribute at next-to-next-to- leading-order). To benchmark the approximation scheme, 
we study electron scattering in both quasi-elastic and pion production channels. 

There have been several experiments measuring the weak response of nuclei across the 



30h37||. which have 



quasi-elastic region to the A excitation peak. In most experiments 
12 C and 16 as the primary target nuclei, the mean energy of the beam is around 1 GeV. 

n 

As emphasized in [2|, we expect our theory to work up to 0.5 GeV, so we do not rely 
on these experiments to constrain the theory at this stage. On the theoretical side, much 



work has been done (e.g., in [29|, l38M60| ) . Most of these papers are based on the global 



1 Recently, a unified framework for handling both coherent and incoherent production has been proposed 
in 



29 



2 In an EFT, there are an infinite number of interaction terms allowed by various constraints. To organize 
them, we can associate power-counting to each vertex and diagram. The calculation can be done in a 
perturbative way by summing diagrams up to some particular power v. See Refs. 

Q A fl Bid for 

detailed discussions about power-counting in QHD EFT. 
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or local Fermi gas approximation and include contributions from one-body currents, with 
improved treatment for final-state interaction (FSI) and A dynamics in the medium. The 
same approach has also been applied in electron scattering (e.g., in 6l|). In 55N60|. scaling 
approaches are used to address quasi-elastic scattering. Moreover, the contribution from 
two-body currents was studied nonrelativistically, for example, in |62(]. In most of these 
calculations, the A dynamics in nuclei is based on the work of [63J , in which the A self-energy 
has been studied using a nonrelativistic model. Parallel to the nonrelativistic studies, some 
work has been initiated in the relativistic framework, QHD EFT, using the local Fermi gas 



(LFG) approximation and inc 
investigated relativistically in 



uding 



69 



one-body currents 



64 



681 ] . The two-body current was 



70j. These works mainly focus on electron scattering. But 



the handling of the A resonance in these papers is somewhat phenomenological. Moreover, 
in both nonrelativistic and relativistic studies, photon production is rarely investigated. 

In this paper, we also apply the LFG approximation [64J to study the one-body current 
contribution. As shown in [ll, |2[, we make use of meson dominance to generate form factors 
for various currents. Because of the built in symmetries in the Lagrangian, conservation of 
vector current and the partial conservation of axial current are satisfied. These properties 
are well preserved in the LFG approximation. Especially for photon production, vector 
current conservation is crucial. The A dynamics, as a key component in this work, is 
explored to some extent. We introduce interactions between A and non-Goldstone meson 
fields to generate the spin-orbit (S-L) coupling that has been introduced in phenomenological 



models 



71 



721 ] . On the other hand, phenomenological knowledge about S-L coupling puts 



constraints on these couplings. Moreover, the A decay width increases in the nucleus, 
because more decay channels are opened up and this effect overcomes the reduction of pion 
decay phase space. Here we follow the phenomenological studies and separate the width 
to the pion decay width and anything else parametrized by the imaginary part of the A 
spreading potential. As a result of opening new decay channels, the flux having excited a 
A resonance can be transferred to channels that do not involve pion or photon production. 
Moreover, Pauli blocking can reduce the pion and photon production cross section further, 
because of the reduction of the final particle's phase space. 3 In this paper, we explore 



3 The binding effect should be important when the neutrino energy is close to threshold, where the simple 
approximations used here are not feasible. But this is clearly not important around 0.5 GeV. 



3 



how both A and nonresonant contributions are reduced compared to those in free nucleon 
scattering. However, we do not include FSI effects for pions and knocked out nucleons. The 
simple treatment can be found in (3, 3' while the complete treatment is implemented 
in various event generators of experiments (e.g., NUANCE l75|), and the GiBUU model 

n □ 

[46|, |47j. Hence we only compare our predictions with the output of NUANCE without FSI. 



The paper is organized as follows. In Sec. [TTJ we first discuss the LFG approximation 
and then apply it to electron quasi-elastic scattering, which serves as a benchmark. In 
Sec. Jill the calculation scheme for pion (photon) production is briefly introduced. Then 
the A dynamics is studied with emphasis on the connection between A-meson interactions 
and S-L coupling. The modification of the A width is also discussed. After that, electron 
scattering at the A peak is studied, and results are compared with data with explanation of 
the missing strength. The cross sections of neutrinoproduction of pions are also shown and 
compared to NUANCE's output. Sec. [TV] is dedicated to NC photon production. Finally 
Sec. [V] contains a short summary. In the appendices, we show detailed kinematic analyses 
for both quasi-elastic scattering and pion production. 



II. QUASI-ELASTIC SCATTERING IN THE LFG APPROXIMATION 



This section serves as an illustration of the LFG approximation used for quasi-elastic 
scattering and for photon and pion production. (See Ref. [2j] for discussion on the free 
nucleon interaction amplitude in all these processes.) Here we make use of the mean- field 
approximation to calculate the nuclear ground state. The relevant leading order Lagrangian 
is 



C = N 



N 



(1) 



^ (dft + igpPv + ig v VA - M + g, 

(where the full Lagrangian can be found in jl, 5| for example). The mean-field approximation 
is presented simply as follows. Inside nuclear matter, vector p 3fl and V 1 , and scalar <ft fields 



4 The predictions from NUANCE shown throughout this paper are obtained from the NUANCE v3 event 
generator [75| . Multiple resonances are considered in NUANCE, but the A dominates. The axial mass 
M\ = 1.10 ± 0.27 GeV is used which is the same as that used by the MiniBooNE experiment for their 



baseline calculations 



data in a separate exercise. 



761 ] . However, the actual backgrounds used in their final analyses were scaled to 
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0.12 




FIG. 1: (Color online) Proton and neutron density in C with Gl and G2 parameter sets. 

develop nonzero expectation values. In the laboratory frame of the matter, only two fields 
{(f) and V°) have nonzero values (but in the isospin asymmetric case, p° can also develop a 
nonzero value). As a result, the nucleon's mass is modified: M* = M — g a (<f>) . At the lowest 
order, the spectrum of nucleons is E(p) = ^/p 2 + M* 2 + g v (V°). Inside a finite nucleus, 
due to different boundary conditions, the mean-field expectation value is space dependent 
and can be calculated numerically. By using this approximation, we can calculate the bulk 
properties of the nucleus, the details of which can be found in Ref. jlj] for example. 

Following [5], we calculate the local density p p / n {f) and field expectation value in 12 C (the 
major nucleus in the MiniBooNE's detector). Figs. [T] and [2] show the results based on Gl 
and Gl parameter sets in We will explore the difference due to the two sets in electron 
quasi-elastic scattering. 

To calculate the electroweak response of nuclei, we use the LFG approximation. This 
approach has been applied in [64J to study electron quasi-elastic scattering. First, by as- 
suming the impulse approximation (IA), the interaction happens every time between probe 
and each individual nucleon This only holds when the transferred momentum is high enough 
that the interference between different nucleons is reduced due to the big recoil. Second, the 
response of the nucleus is the incoherent sum of the response of the fermion gas in different 
regions. This works when the probe's wave length is small enough compared to a charac- 
teristic length scale of the nucleus density profile. The discussion can be summarized in the 
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FIG. 3: The kinematics in the laboratory frame of the nucleus for a lepton interacting with one 
nucleon inside the nucleus. 

following equation: 

In this equation, p ti and pif are the incoming and outgoing lepton momenta, respectively, 
q = pn — pif is the momentum transfer, p n i S{ and p n f Sf are the scattered nucleon's initial 
and final momenta and spin projection, and is the one-body interaction amplitude. The 
kinematic configuration is shown in Fig. [3j The integration over initial and final nucleon 
momenta depends on the space dependent Fermi momentum. A detailed discussion about 
this equation can be found in Appendix [A] 

The interaction amplitude Mfi in Eq. ([2]) can be expressed in terms of various current 

6 



matrix elements (where V?, A* , and Jb m are, respectively, the isovector vector current, 
the isovector axial- vector current and the baryon current, i = ±1,0 

HQ). 

For electron 

quasi-elastic scattering, 

2 

M _ e / jr(Jep) \ / y(/iaci) M\ 
* ~~ EM iil\°EM I ' 

(J^fO = (iV,S|^ + \j%\N,A) , (3) 

where A and S in the state are nucleon isospin. For charged current (CC) quasi-elastic 
scattering, 

M fi = 4v / 2G F K,(jgT ) )(jf a ^) , 
(jCM)*) = ( iV) 5|I(^ + A v )|iV,A) , (4) 

where i — ±1, is the Fermi constant and Vud is the w and d quark mixing element in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. For NC quasi-elastic scattering, 

M fi = AV2G F (J^l)(J ( N h ch , 
{j{ ha d ), } = {N}B \jp _ sin 2 9w r EM \N, A) , (5) 

where 9 W is the weak mixing angle. The electroweak currents of leptons are well known, 
and (N, B\VHJ^ , A l „ )\N, A) can be found in But we need to include the nucleon spec- 
trum modification to the results of j^j, which is straightforward to complete in the LFG 
approximation. 

A short discussion on FSI is in order here. The picture is the following: the interaction 
channels are opened in the initial interacting vertex, and then these channels would couple 
to each other when particles are traveling through the nucleus. The flux among all the 
initial channels are redistributed due to FSI. This picture is adopted in the GiBUU model 
calculations for example for CC and NC processes j^. From conservation of probability, 
assuming the picture mentioned above is valid, we should expect the sum of these channels in 
the initial vertex to match the inclusive data. Moreover, Coulomb distortion of the electron 
is not included in this calculation. 

In the upper panels of Figs. HI El and El we present differential cross sections da/dq°dQ 
for electron scattering off 12 C at given electron energies and scattering angles. In this 
section we only focus on the so-called quasi-elastic peak at the lower energy region, which 
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FIG. 4: (Color online) Inclusive data for differential cross section of electron scattering off 12 C. The 
incoming electron energy is Ei = 0.62 GeV, and the scattering angle is 9\f = 60°. The kinematics 
is measured in the laboratory frame. The data are from Ref. [77|. Explanations of different plots 
can be found in the text. 



is believed to be dominated by one nucleon knock out. The higher energy peak will be 
discussed in Sec. IIHI In Fig. HJ the electron energy is Ei = 0.63 GeV, and the scattering 
angle is Q\f = 60°. The plots "Gl" and "G2" are the calculations done with Gl and G2 
parameter sets j^J. The difference between the two is small. The data are from Ref. 
The validity of the form factors realized by meson dominance needs to be discussed here. In 
this figure, Q 2 ~ 0.3 GeV 2 and |g| ~ 0.55 GeV at the peak. Below the peak, Q 2 is slightly 
higher than 0.3 GeV 2 , and above the peak, Q 2 ^ 0.3 GeV 2 . As discussed in [2], meson 
dominance works when Q 2 ^ 0.3 GeV 2 , and hence it can be applied here. This is also true 
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FIG. 5: (Color online) Inclusive data for differential cross section of electron scattering off 12 C. The 
incoming electron energy is Ei = 0.68 GeV, and the scattering angle is 6\f = 36°. The kinematics 
is measured in the laboratory frame. The data are from [77[. Explanations of different plots can 
be found in the text. 



for Figs. |S]and|nj In Fig. El the electron energy is = 0.68 GeV, and the scattering angle 
is 9if = 36°. In Fig. EJ the electron energy is E { = 0.73 GeV, and the scattering angle is 



6if = 37.1°. The data are from 7t| and 
the "Gl" and "G2" parameter sets. 



78| . Again, we see only small differences between 
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FIG. 6: (Color online) Inclusive data for differential cross section of electron scattering off 12 C. 
The incoming electron energy is E{ = 0.73 GeV, and the scattering angle is Bit = 37.1°. The data 
are from [78j. Explanations of different plots can be found in the text. 

III. PION PRODUCTION 

A. Approximation scheme and A dynamics in the nuclear medium 



By using the LFG approximation detailed before, the formula for the cross section can 
be written as 



a 
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FIG. 7: The kinematics in the laboratory frame of the nucleus for one nucleon inside the nucleus 
interacting with the lepton and producing a pion. 
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FIG. 8: Feynman diagrams for pion production. Here, C stands for various types of currents 
including vector, axial-vector, and baryon currents. Some diagrams may be zero for some specific 



type of current. See Ref. 



2] for the details. 



The details can be found in Appendix [B] The notations for various momenta are explained 
in Fig. [7J All the integrations except the volume integration depend on the space coordinate 
r* through the space dependent Fermi momentum. The amplitude in Eq. fl6]) is similar 
to those in Eqs. ©,01]), and (jSJ), except that the hadronic currents should be changed to 
those relevant to pion production: 

(j( had >) = (N,7r\J ihad ^\N) . 

Several Feynman diagrams contribute here, as shown in Fig. [HI including diagrams with the 
A [(a) and (b)] and all the rest which we define as nonresonant diagrams. See Ref. {2] for 
details about them. Among the medium-modifications of the matrix elements, the behavior 
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of the A needs to be singled out. First, let us focus on the real part of the A self-energy. 
We start from the following Lagrangian (and a similar Lagrangian can be found in js8|): 

b 



-£a;ti-. 



P,V,<t> 



hp / - h v y - m + h s 



A 



bv 



Am Am 



(7) 



Here the A field is given by the Rarita-Schwinger representation, and a,b = ±3/2, ±1/2 
are A isospin indices [l|. At the normal nuclear density, the A is not stable in the nuclear 
medium. So the expectation values of meson fields are not changed in normal nuclei at the 
mean-field level. Similar to the nucleon case, the A spectrum in nuclear matter (without 
the A-pion interaction) is given by 



Pa 



m = 



h v (V°) + ^m* 2 + 0i 
h v (V°)+p*£ 
h v (V°) + ^m* 2 + p% 
m - h s {<j>) . 



The effect of introducing h s and h v couplings on the equation of state (EOS) was analyzed in 



65, 
65, 



80 



811 ] . Some constraints on the couplings, r s = h s /g s and r v = h v /g v , were calculated in 



8jJ. Here we resort to the scattering problem to find other constraints. In pion-nucleus 



scattering studies 
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S-L coupling of the A inside the nucleus was introduced by hand, 
although its origin is not clear in the nonrelativistic model. In this model, a mechanism 
similar to the generation of the nucleon's S-L coupling is used to generate A's. Following 
discussions in |82j and using the Lagrangian in Eq. (J7J), we can estimate the S-L coupling of 
the A as 



2m 2 r dr 
= a(r)S ■ L . 



d (hM+h v (V°})-^±(h v (V )) 

mm r dr v 



S-L 



(8) 



Here, m 



m 



\ (h s {</>) + h v (V )). In Fig. [9j we compare our estimates of a(r) defined in 



Eq. (jo]) with two different phenomenological fits for C. We can see that our estimates based 
on three different parameter sets, r s = 1, r v = 1, 0.9, 0.8, and /„ = —1.0, are consistent 
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FIG. 9: (Color online) The strength of A's S-L coupling in 12 C. Here we compare two phenomenol- 
ogy results with our three calculations based on different parameter sets. The "Horikawa" is from 

i n 

72l | . The "Nakamura" is from [711 ] . All these calculations involve setting the Gl parameter to 
describe the nucleus ground state. We change r v to 1, 0.9, 0.8 while keeping r s = 1, f v = —1. 



with the "Nakamura" result in 7j|, while the "Horikawa" result in 72|] is significantly larger 
than the "Nakamura" result when r > 1 fm. Meanwhile, all the couplings are consistent 
with the "naturalness" assumption, which also motivates our choice of f v = —1. We do not 
show the consequence of r s = r v = 0, since there is no S-L coupling generated in this case. 

Second, we turn to the imaginary part of the self-energy. It is known that Pauli blocking 
effects decrease the width due to reduction of the pion-decay phase space, while the collision 
channels, AiV -h- NN for example, increase the width. The two competing processes have 
been investigated in nonrelativistic models. At normal nuclear density, the net result is to 



increase the width |63|. In phenomenological fits 79] 72|, this increase is taken into account 
by introducing a density- dependent complex spreading potential for the A. Here we follow 
this approach. Above the pion threshold, 



p(r) 



sp 



V x 



P(0) 



(9) 



is the A pion-decay width 



63 



66(. 5 r sp (Vo ~ 80 MeV) is the width in other channels, 



5 In the r,r calculation, only Pauli-blocking is considered. Modifications of the real part of the nucleon and 
A self-energies are not included. 
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and it has been fitted in 7l|, |72(. Below the pion threshold (which is useful in photon 
production), 

p{r) 



sp 



Vo X 



P(0) 



(10) 



n 



;he cross channel of the A diagram, we set the width to zero. Moreover, in the literature 



5Qj, |83], the simple increase of the A width by 5T « 40 MeV has been used for pion 



production: 



T A -> 120 + 40 MeV . 



(11) 



This procedure turns out to work qualitatively, as will be shown later. Furthermore, in 61], 
the A self-energy calculated in 63] is used for inclusive electron scattering off nuclei. In 



Sec. IfflBl we wil 
the width in 61, 



compare our results using Eqs. (EJ) and (jTOj) with those using Eq. ( ITT]) and 
63]. 



B. Pion electroproduction 

Here we focus on the region beyond quasi-elastic scattering in the upper panels of Figs. HJ 
15] and [6] It is believed that the second peak mainly comes from the A excitation inside 
the nucleus. In the upper panels of these figures, we provide our pion-production results 
(without FSI) due to six different calculations. We include the full set of Feynman diagrams 
in the first five calculations, and diagrams with the A in s and u channels in the sixth. 
The difference among the first three calculations is the choice of (r s , r v ) parameter sets: 
(r s = 0, r v — 0), (r s — 1, r v — 1), and (r s = 1, r v = 0.8). In these three, the A width 
shown in Eqs. (]9]) and ( ITU]) is applied. In the fourth calculation, we set (r s = 1, r v = 1) 
and apply the constant shift of the A width as shown in Eq. ( [TT]) . The fifth calculation is 



done by using the A self-energy as calculated in 63|, |84j , which is essentially repeating the 



calculations in 61]. The sixth calculation has (r s = 1, r v = 1) and uses the same A width 
as used in the first three. 

First, let us discuss the location of the A-peak along the q° axis. The different choice of 
r s and r v indicates different binding potentials for A. For (0, 0), the real part of the self- 
energy is the same as in vacuum without any binding. For r s = 1, A has the same attractive 
potential as the nucleon. The vector part tuned by r v provides a repulsive potential. So 
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we can see that (1, 0.8) has a deeper binding potential than (1, 1). Hence, (0, 0) is less 
bound than (1, 1) and (1, 1) is less bound than (1, 0.8). In Figs. HI El and El the location 
of the A peak in first three calculations indeed follows this argument. (We can estimate 
the location of the A peak in a global Fermi gas model. 6 ) The fourth calculation with 
(1, 1) and constant width does not give the correct peak position in the three figures: It 
underestimates (overestimates) A's contribution on the left (right) side of the peak, because 
the constant width assumption overestimates (underestimates) A's width on the left (right) 



side. The fifth calculation by using A modification calculated in [63|, |84| gives the correct 
location of the peak. Comparing the second with the sixth calculation, we can see the 
significance of nonresonant contributions (they use the same set of parameters and the A 
width). 

However, the pion production channel could not explain the full strength of the A peak. 
Meanwhile in the "dip" region between the quasi-elastic scattering and the A peak, the 
calculations also miss strength. This indicates we miss other channels from dip region to the 
A peak. Missing strength at the peak position can be qualitatively explained by considering 
the fourth calculation in the upper panels, whose simple treatment of the A width makes 
analysis transparent. According to Eqs. fl9]) and ffTTT) . we estimate 0.04 ^ r sp ^ 0.08 GeV 
in the sense of averaging over 12 C, and hence 0.08 ^ ^ 0.12 GeV. The comparable 
width of other decay channels shows their importance to the inclusive data. Moreover, 
there are contributions from two-body currents without A as an intermediate state. In the 
lower panels of Figs. HI El and El we add up three different channels: quasi-elastic, pion 
production, and two-body current contributions [labeled as meson-exchange-current (MEC) 
in the plots]. The label "(1, 1, T sp = 0)" for pion production applies under the assumption 
that (r s — 1, r v — 1) and that no new channel takes away the flux from the A pion 
production. The MEC-contributions are from [3] [sjj]. Here the total strength matches 
well with the inclusive data. However a detailed study of different channels in the QHD 
EFT framework is needed to address this issue conclusively. 



6 Following 64j, for 12 C, we assume a global Fermi gas, with the nucleon effective mass M* = Q.75M, 
g s ((j>) w 0.235 GeV, g v (V°) ~ 0.75g s ((f)) (in the laboratory frame), and to* = m — r s g s (cf>) sa 
0.995 GeV(r s — 1). Meanwhile in the s channel, the A momentum is p* A — q + p* H + (1 — r v )g v (V°). It is 
easy to calculate the q°-location of the peak by setting p*£ = to* 2 . For Ei — 0.62 GeV and 9 = 60° (see 
Fig. gl), q° = 0.43 GeV if r v = 1 and q° = 0.39 GeV if r v = 0.8. 
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The difference between our calculations and those in 65j where the QHD model is also 



applied should be mentioned here. Ours are strictly based on the field theory, while in [65 1 
the A is introduced by hand (where they were convoluted with the cross section based on 
a "stable" A theory with a Lorentzian weight function). Moreover, we take into account 
the contribution from other diagrams, which are not considered in 65j. The two results are 
different somewhat, but our choice r s = l,r v = 1, 0.8 is consistent with the analysis in 65 ]. 

Moreover, we can see that the differences in cross sections obtained using r s — l,r v — 
1, 0.8 are not significant, which indicates that the total cross sections of neutrinoproduction 
processes are not sensitive to them either. This will be confirmed by the results in Sec. IIII CI 



C. CC and NC pion production 

Fig. [10] shows the total cross section averaged over proton or neutron number for CC pion 
production in (anti)neutrino- 12 C scattering. We also compare our result with NUANCE's 
output without FSI. In each figure, our calculations including different diagrams and using 
different r s and r v are shown. The "only A" calculation only takes into account A diagrams. 
In the others, all the diagrams up to v = 2 are included. Systematically in all the channels, 
our "only A" calculation is close to the NUANCE output. But other diagrams contained 
in 'V = 2" calculations are not negligible in all the channels around the resonance region, 
especially when the s-channel A contribution is suppressed by the small Clebsch-Gordan 
coefficients (for example, v^ + n — > pT +p + tt°). In the very low energy region away from 
the resonance, nonresonant diagrams dominate. See j^j for the power counting of diagrams. 
Moreover, we check that the contributions of higher order (v ^ 3) diagrams are tiny. We 
also can see that below 0.5 GeV, the (r s = 1 , r v = 0.8) results are bigger than the (1, 1) 
results and the (1, 1) results are bigger than the (0, 0) results. Following the discussions 
in Sec. IIII Bl about the location of the A peak in pion electroproduction, we expect that, at 
a given energy, A excitation occupies more phase space in (1, 0.8) than in (1, 1) and more 
in (1, 1) than in (0, 0). So the pattern among the three different calculations is consistent 
with the qualitative analysis. Here (0, 0) is presented simply for the purpose of comparison, 
and its conflict with the A S-L coupling is presented in Sec. IIII Al 

One question needs to be raised: Do we have A-dominance in the nuclear scattering 
around 0.5 GeV ? If we compare the 'V = 2" calculations with "only A" calculations in every 
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FIG. 10: (Color online) Total cross section per proton or neutron for the incoherent CC pion 
production in neutrino- and antineutrino- 12 C scattering. 



channel, the answer is no. It turns out that the A contribution is strongly reduced due to the 
broadening of its width, compared to its contribution in free nucleon scattering. Meanwhile 
nonresonant contributions are reduced by Pauli blocking. To see this qualitatively, compare 
our results here with the cross sections shown in yf for production from free nucleons. In 

n 

[2[, two different calculations can be found, including "Only A" and ' v = 2". We just show 
the total cross sections at E v = 0.5 GeV in Tab. |T] for neutrino scattering. For example, 
"p,p7r + " indicates the channel u+p — > +p + 7r + . "(f)" and "(b)" correspond to scattering 
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<j(Kr 39 cm 2 ) 


only A (f) 


v = 2 (f) 


Nonresonant (f) 


Only A (b) 


v = 2 (b) 


Nonresonant (b) 


p,pir + 


0.56 


0.85 


0.29 


0.33 


0.48 


0.15 


n, mr + 


0.088 


0.105 


0.017 


0.056 


0.060 


0.004 


n,pir° 


0.117 


0.258 


0.141 


0.069 


0.153 


0.084 



TABLE I: Total cross sections averaged over number of protons or nucleons for CC pion production 
in neutrino- 12 C scattering at E u = 0.5 GeV. See the text for detailed explanations. In the nuclear 
scattering, r s = r v = 1. 



from free nucleons and from bound nucleons in 12 C respectively. In both "only A (b)" and 
u u = 2 (b)," r s = r v = 1 (and note that calculations with only A and r s = r v = 1 are 
not shown in the figures). "Nonresonant (b)" is the difference between the two, and can be 
viewed qualitatively as the contributions of the nonresonant diagrams. 7 The labeling for 
free nucleon scattering is the same. We can see that the A contribution in nuclear scattering 
has been reduced systematically by around 50% in all channels, compared to its contribution 
in nucleon scattering; the nonresonant contributions are also strongly reduced. Clearly, the 
nonresonant contributions are not negligible in both nucleon and nuclei scattering. The 
same situation occurs in the antineutrino scattering channels and hence are not shown 
explicitly. This underscores the importance of including nonresonant contributions in CC 
pion production. 

In Fig. [TlJ we show the total cross section for NC pion production from 12 C. The 
categorization of the different calculations are the same as those for CC scattering. Again 
the NUANCE output is close to our "only A" calculation. Among the first three calculations 
in each channel, at fixed (anti)neutrino energy, (1, 0.8) gives a larger cross section than (1, 1) 
and (1, 1) gives a larger cross section than (0, 0). This is the same as in the CC production, 
which has been explained in terms of kinematics. Moreover, we can see how A-dominance 
is violated in the NC CclS6 ; clS shown in Tab. [TT] (in which the labelings are the same as 
those in Tab. [H and free nucleon scattering results are from {2]). The same is true for 
antineutrino-nucleus scattering. 



7 In principle, there are interferences between contributions from A and other diagrams. At E w = 0.5 GeV, 
we can assume in most of phase space that the A is "on shell" while contributions from other diagrams 
are real, and hence the interferences are small. 



18 



0.25 
0.2 

0.15 
0.1 

0.05 




only A, (0, 0) ' 


I I 


v=2, (0, 0) 




v=2, (1, 1) 




- v=2, (1, 0.8) - 




NUANCE 






V p->V p 71 




-tccrrfW*'"" 'i |_ 





0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

E v (GeV) 



E 

o 

CO 

I 

o 



(b) 0.3 
0.25 
0.2 
0.15 
0.1 
0.05 




only A. (o. 0)' 

v=2, (0, 0) 

v=2, (1, 1) 

v=2, (1, 0.8) 

NUANCE ■■ ■■ 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

E v (GeV) 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 
E v (GeV) 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

E v (GeV) 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

E v (GeV) 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

E v (GeV) 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 

E v (GeV) 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 



E v (GeV) 



FIG. 11: (color online). Total cross section per proton or neutron for the NC pion production in 
neutrino- and antineutrino- 12 C scattering. 
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<j(Kr 39 cm 2 ) 


only A (f) 


v = 2 (f) 


Nonresonant (f) 


Only A (b) 


v = 2 (b) 


Nonresonant (b) 


p,pir° 


0.194 


0.230 


0.036 


0.101 


0.121 


0.020 


n, nir° 


0.194 


0.234 


0.040 


0.101 


0.123 


0.022 


n,p-ir~ 


0.089 


0.149 


0.060 


0.045 


0.082 


0.037 


p, nir + 


0.089 


0.155 


0.066 


0.045 


0.088 


0.043 



TABLE II: Total cross sections averaged over number of proton or nucleon for NC pion production 
in neutrino- 12 C scattering at E v = 0.5 GeV. See the text for detailed explanations. In the nuclear 
scattering, r s = r v = 1. 

IV. NC PHOTON PRODUCTION 



In this section, we study NC photon production from 12 C. The calculation is done in the 
same way as in pion production, except that the hadronic current in Eq. ([5]) is changed to 
the following: 

(jW") = (N,-f\J ihad)ll \N) . 

The Feynman diagrams are the same as those in Fig. [8] with the final ir line substituted by 
the final 7 line. See Ref. [2] for detailed discussion about them. Again we need to implement 
the change of the baryon spectrum when we apply the formula in as we do in previous 
calculations. Because of built in symmetries in our model, conservation of the vector current 
is automatically satisfied, which is important for photon production. The difference in the 
kinematic analysis, compared to that in pion production, is due to the zero mass of the 
photon. Moreover, we apply an energy cut on the photon energy in the laboratory frame, 
Ej ^ 0.15 GeV, motivated by the MiniBooNE's detector efficiency. This also eliminates the 
infrared singularity and simplifies the calculation. 

In Fig. [12], the total cross sections averaged over proton or neutron number are shown. 
Four different calculations are compared. The first "only A" is the same as before. u u = 3" 



2 contact diagrams 

contribute, and there are only two v = 3 contact vertices contributing (See Ref. 2J for 



calculations include all the v ^ 3 diagrams. It turns out no v 
contribu 
details): 



As we have checked, the contributions of these two are small compared to those of the A and 
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FIG. 12: (Color online) Total cross section per proton or neutron for the NC photon production 
in the neutrino- and antineutrino- 12 C scatterings. Our calculation is done with a photon energy 



cut EL > 0.15 GeV. 



existing nonresonant diagrams, which should be expected according to the power-counting. 
Here, we have assumed their strength are due to both the u and p meson anomalous inter- 
ac tl o n ve rtlCeS ( Cl = 1.5 and „ = 0.8) fi Q. Moreover tor these captions, e llmglng 
r s and r v does not change the total cross section significantly, which is also observed in 
the differential cross section for pion electroproduction. In the three v = 3 calculations for 
different channels, (1, 0.8) gives a bigger cross section than (1, 1) and (1, 1) is bigger than 
(0, 0). This pattern has been explained in pion production. We also see that the NUANCE 
output is close to the "only A" calculation and smaller than the full calculations, which 
should be expected from the comparison in pion production. 

In addition, in Tab. Ill II we show how the A significance changes from neutrino-nucleon 
scattering to neutrino-nuclus scattering (free nucleon scattering results are from [2| with a 
change on photon energy cut: E 1 ^ 0.15 GeV): the A contribution is strongly reduced, and 
the nonresonant contribution is reduced less significantly. Since we have put a constraint on 
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<j(l(r 42 cm 2 ) 


only A (f) 


v = 1 (f) 


Nonresonant (f) 


Only A (b) 


v = 1 (b) 


Nonresonant (b) 




1.89 


2.49 


0.60 


0.98 


1.50 


0.52 


n, 717 


1.89 


2.25 


0.36 


0.97 


1.24 


0.24 



TABLE III: Total cross sections averaged over number of proton or nucleon for NC photon produc- 
tion in neutrino- 12 C scattering at E u = 0.5 GeV. Here ^ 0.15 GeV for both types of scattering. 
In the nuclear scattering, r s = r v =l. 



the minimum photon energy, the lower energy events are not included in the results and the 
Pauli blocking effect is not significant. That explains why the nonresonant contribution is not 
quite suppressed. And the reduction of the A contribution is mainly due to the broadening 
of its width. We also expect the Pauli blocking effect to be less significant with higher energy 
neutrinos. Furthermore, the same pattern about the reduction of cross sections happens in 
antineutrino scattering. Based on Tab. IIII[ we need to include nonresonant contributions in 
photon production, as emphasized in pion production. 



V. SUMMARY 



Neutrinoproduction of photons and pions from nuclei provides an important background 
in neutrino-oscillation experiment and must be understood quantitatively. Especially, we 
are interested in the possible role of NC photon production in the excess events seen in the 
MiniBooNE experiment at low reconstructed neutrino energy. In Ref. [2|, we have calibrated 
our theory — QHD EFT with A introduced — by calculating photon and pion production from 
free nucleons up to E u = 0.5 GeV. In this work, the theory is applied to study the produc- 
tion from nuclei. Here we make use of the LFG approximation and Impulse Approximation, 
and include only one-body current contributions. In the mean-field approximation of the 
nuclear ground state, the change of the baryon spectrum is represented by introducing an 
effective mass for baryons, which leads to the change of one-body currents in this calculation. 
The calculation for electron quasi-elastic scattering and electroproduction of pion serves as 
a benchmark for our approximation schemes. We then proceed to calculate the neutrino- 
production of pion and photon from 12 C, and show the plots for total cross section in every 
channel. First, we present calculations for pion production up to next-to-leading-order with 
different r s and r v parameters as constrained by the phenomenological study. It turns out 
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that total cross sections are not very sensitive to changes of these parameters. Then in NC 
production of photon, although we show the result up to v = 3 order, there are no v = 2 
contributions from contact terms, and as we have checked already the v = 3 contributions 
due to ci and ei, related to the so-called anomalous interactions, are tiny (the same has been 
shown for nucleon scattering in |2j). Again, the total cross section of photon production is 
not sensitive to choice of different r s and r v . In all the plots, the A contributions are singled 
out and compared with the full calculations. Moreover, we also compare our results with 
the output from NUANCE, and we find that the NUANCE output is close to our "only A" 
calculation with (r s — 0, r v = 0) for both pion and photon production, which should be 
expected since the A dominates in NUANCE. 

In the calculation, the A dynamics in nuclei is a key component. The dynamics has been 
investigated in a nonrelativistic framework and also initiated in the QHD model. Parallel 
to the modification of the nucleon's spectrum, the A-meson couplings (related to r s and 
r v ) introduced in our theory dictates the real part of the A self-energy. The couplings 
are used to explain the S-L coupling of A. Meanwhile the phenomenological result about 
S-L coupling based on nonrelativsic isobar-hole models puts an interesting constraint on 
the A-meson coupling strengths, which is complementary to the constraints based on an 
EOS consideration. The A width is treated in a simplified way, as we take advantage of 
the existing result that shows an increase of the width due to the opening of other decay 
channels. In pion electroproduction, the pion-production (without FSI) result gives a correct 
prediction for the location of the A-peak. We argue that this deficit is due to the absence 
of other channels. By adding contributions from two-body currents (from other relativistic 
studies) to our quasi-elastic and pion production (and turning off A broadening), we can 
explain the inclusive electron scattering strength. The investigation on A dynamics and 
two-body currents, which plays an important role in nuclear response and other problems, 
certainly needs to be pursued further in QHD EFT. 

Moreover, because of the broadening of the A width, we expect that in both pion and 
photon productions, the A contribution is much less in nuclear scattering than in nucleon 
scattering. But the reduction of nonresonant contributions would be less at higher energies 
(beyond 0.5 GeV), because the Pauli blocking effect should be less important. In Tabs. HI HT1 
and IHIl we have shown explicitly the cross sections at E v = 0.5 GeV due to A and nonreso- 
nant contributions in both neutrino-nucleon and neutrino-nucleus scattering. Although we 
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see the reduction of nonresonant contributions for pion production in Tabs. |T] and [TTJ, we see 
a smaller reduction for photon production in Tab. IIII1 This is consistent with the picture 
that the nonresonant contribution is reduced because of Pauli blocking. The same situa- 
tion occurs in antineutrino scattering. This conclusion is important for future investigations 
of higher energy neutrino scattering, which may be relevant to MiniBooNE's excess event 
problem. 

Since our calculation is based on a QHD EFT Lagrangian with all the relevant symmetries 
built in, conservation of vector current is manifest. This is crucial for photon production. 
Also partial conservation of the axial current is a necessary constraint in the problem. By 
using the mean-field approximation and the LFG model, these constraints are satisfied in a 
transparent way. 

We are currently working on coherent pion and photon production from nuclei by applying 
this QHD EFT, which may also be relevant to the MiniBooNE low energy excess event 
problem. 
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Appendix A: kinematics for quasi-elastic scattering 

The analysis of the kinematics is for scattering from nuclear matter, and can be easily 
generalized in the LFG model. The kinematic variables are shown in Fig. [5J and discussed 
following Eq. <^I§. From the mean- field theory in QHD EFT, we know that the leading 
order Hamiltonian gives rise to the nucleon spectrum in nuclear matter as p° = g v (V°) + 
^M* 2 + p2, M* = M - g s (4>). Then we can define p* n ° = p° - g v (V°) = V 'M* 2 + p* = 
Vm* 2 + (f n y. This can be generalized from the laboratory frame to an arbitrary frame. In 
the LFG model, we consider each neighborhood inside the nucleus as a homogeneous system; 
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the field expectations, (<j>{x)) and (V^^x)), are space-time dependent (and in the laboratory 
frame, they only depend on the space coordinate). In the following, we always work in the 
nuclear laboratory frame. The covariance of our calculation is more transparent with the 
p* n ^ variables than with p%. For example energy momentum conservation is q + p* ni = p* x. 

Next we derive the formula for the total cross section. Suppose Mfi is the covariant 
interaction amplitude between the probe and each individual nucleon with specific initial 
and final states. We have 



a 



Pauli blocking leads to constraints on the integration of p* ni and i.e. \p* ni \ ^ pp and 
\p^,f\ ^ Pf- Here pp is the Fermi momentum related with the local density. The two 
constraints can be expressed by using factors 0[pp+p* n 2 ~ (.Pni ' V) 2 /V 2 } an d 0[—p F —p* n 2 + 
(Pnf ' V) 2 /V 2 ]- in the following, we will not include them explicitly. We know that 



d 3 pXi d 3 f?nf r-4/ , * * \ f 11 i *0 1 

2p *0 2p*° f W+Pni-Pnf) ~ J d <Pp* n i d Pni^ I <• >s. ^//-;„)-( 2,,' >„','+./- I r>|,7l|/7;„ i ' 



By using this, we have the total cross section as 



a 



, 1 f d 3 p lf dp^d<j> KiSr 



,0. 



Meanwhile to make our phase space analysis simple, we can integrate over d\q\ and dq 

° = j ^ d ^JPni d ^ l4^l^' M/lf • (A3) 

Sf,Si 

Now, we need to calculate the boundary of the phase space in Eq. (1A3I) . From the lepton 
kinematics, we can determine the boundary of \q\ (p^ = 



1. = \pu\ + \lEl-Mf fl (A4) 



l$U = Wu\ ~ ^El - MJ . (A5) 
For a given we have the following constraints based on the lepton kinematics: 



q° < E H - (E lf ) mm = E u - ^{\p u \ - |gl) 2 + Mf f , (A6) 
q° > E u - {Eif) max = . (A7) 
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However, there are further constraints on q° for a given \q\ due to the hadron kinematics. 
For a given set of q°, \q\, cos(Zgj3*j) = (2g°p*° + g 2 )/(2|<f||p^|) has to be physical. This 
requires 



cos(lqp* ni )\ < 1 



M 

2 



4M 



(A8) 



Eq. (IA8|) gives a lower bound of which is also required to be below the Fermi surface: 
\Prii\ ^ Vf- Combining pp ^ p~ and the constraints in Eqs. flA6j) and (jA7|), we find 



max 



mm 



\q\-p F ) 2 + M* 2 -E F ,0 



+ p F f + M* 2 - Ep-i Eu - x /(\p u \ - \q]) 2 + Mff 



(A9) 



(A10) 



Moreover, the constraint \p^f\ ^ Pf is not present in the former discussion, but is taken care 
of in the numerical calculation. 



Appendix B: kinematics for pion production 

The kinematic variables are defined in Fig. [7] in the laboratory frame. Except for the 
7r momentum k n , all the others are defined in Appendix [K\ The variables defined in other 
frames will be mentioned explicitly. First we have 

( dV — I ^ ^ d ^ lf ^ 
a J 2pl J (2tt)>;° (27r)32A;0 (2^2^ (2k)*2j% 

x(27r)V(g +p* ni -p* nf - K)J2\M fl \ 2 • 

Sf,Si 

The constraints on p* ni and p* n f, i. e. |p^| ^ p F and |p^J ^ Pf, are always implicit in the 
formula. 

One way to think about the phase space as follows: Given specific values for q and p* nil 
the final pion and nucleon invariant mass M nn are fixed, and then the degrees of freedom 
in the isobaric frame (final pion and nucleon's center-of-mass frame) is the angle of kj n , i.e. 
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fir . So we have 

K Iir 



(27r) 3 2p;° / (2tt) 3 2££ 
1 |fc /7r 



(27r)V(g + K,-K/-Av) 

*[(? + j0 2 -<J . 



= / dM nn dflt 

J fc ^(2vr) 2 2 

In the above, we have made use of the following identities: 



dM m 



M* I + \kj^ + sjM^ + \k I ^ , 

2 1 I ■ '- i 1 /2 



.1/- - .1/*- + M 



2M» 



-^/n/ — M nn — E l7r 



E 



Inf 



M, 



Then analogous to the analysis in the quasi-elastic scattering case, we have 



Pni 



(2vr) 3 2^ 



S[(q + P* ni T - M, 



I 2 1 

wni 



-s 



3 lcos(Z^)=(2gO p; o +g 2 +M ,2_ M 2 n)/2k1K j . 



So finally: 



/ dvdJ>dWM„d P '°d<i> F jn k 



(2>r) 7 32(p«) 



Next, we need to determine the boundary of phase space in terms of these variables. 
First, it is clear that no constraint needs to be applied to fi^ and 4>jr ni - Second, for a given 
set of q,q°, and M wn , to make sure | cos(Aqp* nj )\ < 1, there is a constraint on besides 

-M\f m \ -q 2 < + m* 2 - Ml <: 2\q\\fj - q 2 , 



m\ + ^W 2 + ^ - { -^(q°) 2 > , 



x = Mj-M' 2 



Ml„-M* 2 



\fni\ >P = 



(Bl) 
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Third, for a given set of q, q°, there is a constraint on M nn , such that p ^ p F . From 
M 2 n = q 2 + M* 2 + 2q°p* n ° l - 2\q\\fJ cos(Zqp* ni ), we have 

q 2 + M* 2 + 2q°E F - 2\q\p F <: M 2 n «C q 2 + M* 2 + 2q°E F + 2\q\p F . 

And to open the pion production threshold, we need 

+ M*) 2 ^ M 2 n , 

So, we have 

max((M w + M*) 2 , g 2 + M* 2 + 2q°E F - 2\q\p F ) 

< M 2 n ^ g 2 + M* 2 + 2q°E F + 2\q\p F . (B2) 

Fourth, for a given |g|, there is a constraint on q° such that (M nn ) min ^ (M 7rn ) maa ,. We have 

g 2 + M* 2 + 2q°E F + 2\q\p F ^ (Af w + M*) 2 



^ q° ^ max(y/(M„ + M*) 2 + (\q\~ Pf) 2 - E F , 0) . (B3) 

However, there are further constraints on q° due to the lepton kinematics, which has been 
shown in Eqs. (1A6[) and (1A7|) . Together with Eq. (1B3I) . we have the boundary of g°: 



q° max = E H -J(\p u \-\q\) 2 + M 2 f , (B4) 



q° ■ 

lrnm 



max( v/(M w + M*) 2 + (|g1 - p F ) 2 - E F , 0) . (B5) 



Eqs. ( 1A4I) and ( 1A5I) give constraints on |g|. And there are further constraints due to hadron 
kinematics. We have to make sure that q Q max ^ Qmin- But it is complicated to obtain an 
analytic expression for |g| based on this constraint. In the numerical calculations, we made 
use of another boundary by assuming a static nucleon in vacuum. Then, we can say the 
allowed region of |<f| is always inside the previous region. Solving g^ aa , ^ q Q min with p F = 
and M* = M gives us: 

lol < - Eu+M (B6) 

m ^ 2(1 ' W 

® > — w^m • (B7) 
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In the above, 

R - Eh 

~ e^Tm 



M A = ^{Eu + Mf-El, 

\M\ + (Mn + Mf - Mf f f 
A = [ \e k + m)2 — ~ ~ (1 2 a)(M„ + M)< 



So, Eqs. (JA4J), (JA5J), (IB6|) and (JB7J) are the bounds used in the numerical calculations. And 
to map out the physical region, we simply try and check. It is also complicated to determine 
an analytic expression for the threshold value of En in the LFG model. However, it is simpler 
to work out the value for pion production off a static nucleon, which is 



En > 



(M n + M + Miff - M 2 
2M 



So, the difference between the true threshold and the value calculated above is essentially 
the binding energy. 
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